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Abstract 

We review the motivations for the search of stochastic backgrounds of gravitational 
waves and we compare the experimental sensitivities that can be reached in the near 
future with the existing bounds and with the theoretical predictions. 
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1 Motivations 



A possible target of gravitational wave (GW) experiments is given by stochastic backgrounds 
of cosmological origin. In a sense, these are the gravitational analog of the 2.7 K microwave 
photon background and, apart from their obvious intrinsic interest, they would carry extraor- 
dinary informations on the state of the very early universe and on physics at correspondingly 
high energies. To understand this point, one should have in mind the following basic physical 
principle: 

• a background of relic particles gives a snapshot of the state of the universe at the time 
when these particles decoupled from the primordial plasma. 

The smaller is the cross section of a particle, the earlier it decouples. Therefore particles 
with only gravitational interactions, like gravitons and possibly other fields predicted by 
string theory, decouple much earlier than particles which have also electroweak or strong 
interactions. The condition for decoupling is that the interaction rate of the process that 
mantains equilibrium, F, becomes smaller than the characteristic time scale, which is given 
by the Hubble parameter H, 



(we set ^ = c = 1). A simple back-of-the-envelope computation shows that, for gravitons. 



so that gravitons are decoupled below the Planck scale Mpi ~ 10^^ GeV, i.e., already 10"'^'^ 
sec after the big-bang. This means that a background of GWs produced in the very early uni- 
verse encodes still today, in its frequency spectrum, all the informations about the conditions 
in which it was created. 

For comparison, the photons that we observe in the CMBR decoupled when the temper- 
ature was of order T ~ 0.2 eV, or 3 ■ 10^ yr after the big bang. This difference in scales 
simply reflects the difference in the strength of the gravitational and electromagnetic inter- 
actions. Therefore, the photons of the CMBR give us a snapshot of the state of the universe 
at t ~ 3 ■ 10^ yr. Of course, from this snapshot we can also understand many things about 
much earlier epochs. For instance, the density fluctuations present at this epoch have been 
produced much earlier, and the recent Boomerang data indicate that they are compatible 
with the prediction from inflation. Therefore, from the photon microwave background we 
can extract informations on the state of the Universe at much earlier times than the epoch of 
photon decoupling. In this case, conceptually the situation is similar with trying to under- 
stand the aspect that a person had as a child from a picture taken when he was much older. 
Certainly many features can be inferred from such a picture. Gravitational waves, however, 
provide directly a picture of the child, and therefore give us really unique informations. 

It is clear that the reason why GWs are potentially so interesting is due to their very 
small cross section, and the very same reason is at the basis of the difficulty of the detection. 
In the next sections we will review the experimental aspects of the search for stochastic 
backgrounds and the theoretical expectations. Details and more complete references can be 
found in ref. 



T<^H 



=^ decoupled 
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2 Characterization of stochastic backgrounds of GWs 

The intensity of a stochastic background of GWs can be characterized by the dimensionless 
quantity 

^^if) - - , (3) 

Pc dlogf 

where pg„ is the energy density of the stochastic background of gravitational waves, / is the 
frequency {u = 27if ) and pc is the present value of the critical energy density for closing the 
universe. In terms of the present value of the Hubble constant Hq and Newton's constant 
Gat, the critical density is given by Pc ~ 3ifg/(87rGjv)- The value of Hq is usually written as 
Hq = ho X 100 km/(sec-Mpc), where parametrizes the existing experimental uncertainty. 
However, it is not very convenient to normalize Pg^ to a quantity, Pc, which is uncertain: 
this uncertainty would appear in all the subsequent formulas, although it has nothing to 
do with the uncertainties on the GW background itself. Therefore, we rather characterize 
the stochastic GW background with the quantity /iQfig„(/), which is independent of ho- All 
theoretical computations of a relic GW spectrum are actually computations of dp^^ / dlog f 
and are independent of the uncertainty on Hq. Therefore the result of these computations 
is expressed in terms of /iQ^^gw, rather than of Qgw. 

To understand the effect of the stochastic background on a detector, we need however 
to think in terms of amplitudes of GWs. A stochastic GW at a given point x — Q can be 
expanded, in the transverse traceless gauge, as 

M^)= E r df [dQhA(f,n)e-'-'f'ei{Q), (4) 

A=+,x ^ 

where hA{—f,^) = h\{f,Cl). Q is a unit vector representing the direction of propagation 
of the wave and dCl — d cos 9d(f). The polarization tensors can be written as e'^{d) — 
TT^a'iTT'b ~ nafi}, , e^^(Q) = rhafib + TiaThb , with m, fi unit vectors ortogonal to Cl and to each 
other. With these definitions, e^i^{Vt)e^' '""^ [Vt) = 25^"^' . For a stochastic background, assumed 
to be isotropic, unpolarized and stationary, we can define the spectral density Sh{f) from the 
ensemble average of the Fourier amplitudes, 

{h*A{f, ^)hA'if', ^')) = sif - f)^s\n, n')5AA'lsu{f) , (5) 

where S'^{Q,fl') = S{(f) — (j)')S {cos 9 — cos 9') . Sh{f) has dimensions Hz"-*- and satisfies Sh{f) — 
Sh{—f)- The choice of normalization is such that 



/ 



dndn' {h*A{f, n)hA'{f\ n')) = 6{f - f')6AA'lSh{f) . (6) 



figw and Sh are related by 



4-^ .3 



^gw(/) = ^/Skif) . (7) 

In general, theoretical predictions are expressed more naturally in terms of hlflgviif), while 
the equations involving the signal-to-noise ratio and other issues related to the detection are 
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much more transparent when written in terms of Sh{f)- Eq. (J^) is the basic formula for 
moving between the two descriptions. 

The characteristic amphtude hc{f) is instead defined from 

/•/=oo 

(h^tm^'it)) = 2 dOogf) hiu). (8) 

Jf=0 

hc{f) is dimensionless, and represents a characteristic value of the amplitude, per unit log- 
arithmic interval of frequency. The factor of two on the right-hand side of eq. (H) is part of 
the definition, and is motivated by_ the fact that the left-hand side is made up of two contri- 
butions, given by {h*^h^) and {h*^hy). In a unpolarized background these contributions are 
equal, while the mixed term {h\hy) vanishes, eq. (^. The relation betheen he and Sh is 

hlU) = 2fS,{f) . (9) 

Actually, hc{f) is not yet the most useful dimensionless quantity to use for the comparison 
with experiments. In fact, any experiment involves some form of binning over the frequency. 
In a total observation time T, the resolution in frequency is A/ = 1/T, so one does not 
observe hlQg„{f) but rather 

^^^"^^ d{\og f) hln,M) ^ ^Kn^if) , (lo) 

and, since hlQg„{f) ~ f^Sh{f) ~ Ph^f), it is convenient to define 

h,{f,Af) = K{f)(^^^ ' . (11) 

Using 1/(1 yr) ~ 3.17 x 10^^ Hz as a reference value for A/, and 10~® as a reference value 
for /ig^gw, one finds 



A/) = 7.111 X 10- (i^y c^hMiy ( ^ 



3/2 /!,2o 1/2 / A / \ 1/2 

^3.17 X 10-8 Hz ' ^ ' 

3 Characterization of the detectors 

The response of a detector is characterized by two important quantities: the strain sensitivity 
hf, which gives a measure of the noise in the detector, and the pattern functions F^(9,(l)), 
which reflect the geometry of the detector. 

3.1 Strain sensitivity 

The total output of the detector S{t) is in general of the form 



S{t) = s{t)+n{t) 



(13) 
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where n(t) is the noise and s{t) is the contribution to the output due to the gravitational 
waves. Both are taken to be dimensionless quantities. If the noise is gaussian, the ensemble 
average of the Fourier components of the noise, n(/), satisfies 

(r(/K/')) = 5(/-/')^5„(/) (14) 

(actually, non-gaussian noises are potentialy very dangerous in GW experiments, and must 
be carefully minimized or modelled). The above equation defines the function Sn{f), with 
Sni—f) = Snif) and dimensions Hz~^. 

The factor 1/2 is conventionally inserted in the definition so that the total noise power is 
obtained integrating Sn{f) over the physical range < / < oo, rather than from — oo to oo, 

{n\t))= / d/5„(/). (15) 

The function Sn is known as the spectral noise density. This quantity is quadratic in_the 
noise; it is therefore convenient to take its square root and define the strain sensitivity hf 

hj^^/sJT), (16) 
where now f > 0; hf is linear in the noise and has dimensions Hz~^/^. 



3.2 Pattern functions 

GW experiments are designed so that their scalar output s{t) is linear in the GW signal: if 
hab{t) is the metric perturbation in the transverse-traceless gauge, 

s{t) = D'^'Kbit) . (17) 

D^^ is known as the detector tensor. Using eq. (^), we write 

df / dn hAif, n) e-^-'f' D'^'eiin) . (is) 

/i=+,x -°° 

It is convenient to define the detector pattern functions Fa{^), 

Fa{Ci) = D-'ciiQ) , (19) 

SO that 

df / dQhA{f,n)FA{n)e-'-'^\ (20) 
and the Fourier transform of the signal, §(/), is 

5(/)= E I d(lhA{fMFA{(l). (21) 

Explicit expressions for the pattern functions of various detectors can be found in ref. 0]. 
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3.3 Single detectors 

For a stochastic background the average of s{t) vanishes and, if we have only one detector, 
the best we can do is to consider the average of Using eqs. @ and (pO|), 



oo 



{s\t)) = F / dfS.if), (22) 
Jo 

where 

p^jd^ F\h)F\h) (23) 

is a factor that gives a measure of the angular efficiency of the detector. For interferometers 
F = 2/5, while for cilindrical bars F = 8/15. 

In a single detector a stochastic background will manifest itself as an eccess noise. Com- 
paring eqs. ([l5|) and eq. (^2|) , we see that the signal-to-noise ratio at frequency / is 



SNR 



Snif) 



(24) 



We use the convention that the SNR refers to the GW amplitude, rather than to the GW 
energy. Since hlQgv^if) ~ ^c(/) ^^e SNR for the amplitude is the square root of the SNR 
for the energy, hence the square root in eq. (p^. 

Using eq. (|^ and Sn{f) = hj-, we can express this result in terms of the minimum 
detectable value of h^Qg^, at a given level of SNR, as 

h^n-^^( f) - 10- (^] ' ( ^ ' (25) 

V^gwUj-lU ^ I lOOHzi llO-22Hz-i/2 ' ■ ^ ' 



3.4 Correlated detectors 



A much better sensitivity can be obtained correlating two (or more) detectors. In this case 
we write the output Si{t) of the i—th detector as Si{t) = Si{t) + riiit), where i = 1,2 labels 
the detector, and we consider the situation in which the GW signal Sj is much smaller than 
the noise nj. We correlate the two outputs defining 



/.T/2 /.T/2 

Si2= dt dt' Si{t)S2{t')Q{t-t'), (26) 

J-T/2 J-T/2 



where T is the total integration time (e.g. one year) and Q is a real filter function. The 
simplest choice would be Q{t — t') = 6{t — t'). However, if one knows the form of the signal 
that one is looking for, i.e. Sh{f), it is possible to optimize the form of the filter function. 
It turns out that the optimal filter is given by ||, |^, |^, P] 

~ ^X/)WL (27) 

with c an arbitrary normalization constant; Sh{f) is the spectral density of the signal and 
^nKf) the noise spectral densities of the two detectors. Since in general Sh{f) is not known 
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a priori, one should perform the data analysis considering a set of possible filters. For 
many cosmological backgrounds, a simple set of power-like filters Sh{f) ~ should be 
adequate. Actually, many of the most interesting cosmological backgrounds are expected to 
be approximately fiat within the window of existing or planned experiments. 

The function r(/) in eq. (|27|) is the (unnormalized) overlap reduction function; it is 
defined as 



exp 1 2TTifCl • — I , (28) 



c 



. A 

where Ax is the separation between the two detectors. It gives a measure of how well the 
two detectors are correlated; for instance, if the distance between them is much bigger than 
the typical wavelength of the stochastic GWs, the two detectors do not really see the same 
GW backgrounds, and noting is gained correlating them. It also depends on the relative 
orientation between them. 

It is conventional to introduce the (normalized) overlap reduction function ■j{f) [Q, |] 

T(f) 

7(/) = ^, (29) 

where 

^12 = / ^E^i'(^)^2'(^)laligned. (30) 

Here the subscript means that we must compute F12 taking the two detectors to be perfectly 
aligned, rather than with their actual orientation. If the two detectors are of the same type, 
e.g. two interferometers or two cylindrical bars, F12 is the same as the constant F defined 
in eq. (^3]). The use of r(/) is more convenient when we want to write equations that hold 
independently of what detectors (interferometers, bars, or spheres) are used in the correlation 
(furthermore, in the case of the correlation between an interferometer and the scalar mode 
of a sphere F12 = 0, so this normalization is impossible; then, one just uses r(/), which is 
the quantity that enters directly in the SNR). 

The SNR for two correlated detectors, with optimal filtering, is 



SNR 



2T dfT\f) 
Jo 



slU). 



(31) 



Comparing with eq. ( plf ) we see that now the SNR is given by an integral over all frequencies, 
rather then by a comparison of Sh{f) with Sn{f) at fixed /. We denote by A/ the frequency 
range where both detectors are sensitive and at the same time the overlap reduction function 
r(/) is not much smaller than one. Then, if the integration time T is such that TAf ^ 1, 
correlating two detectors is much better than working with the single detectors. 

Finally, it can be useful to define a dimensionless characteristic noise associated to the 
correlated detectors. It turns out that the quantity that is meaningful to compare directly 
to the characteristic dimensionless amplitude hc{f) of the GW (defined in sect. |^) is a 
characteristic noise hn{f) defined by 

hnif) = arlfy/^ (t^t) ' ^^^^ 
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The characteristic noise hn{f) is a useful but approximate concept, since the correct expres- 
sion of the SNR for correlated detectors is obtained integrating over all frequencies, eq. (|3T|) , 
and not comparing and hc{f) at fixed /. 



4 Experiments 
4.1 Resonant bars 

There are five ultracryogenic resonant bars around the world, in operation since a few years: 



detector 


location 


taking data since 


NAUTILUS 


Frascati, Rome 


1993 


EXPLORER 


Cern (Rome group) 


1990 


ALLEGRO 


Louisiana, USA 


1991 


AURIGA 


Padua, Italy 


1997 


NIOBE 


Perth, Australia 


1993 



Table 1: The existing resonant bars. 



Compared to interferometers, they are narrow band detectors, since they operate in a 
frequency band of the order of a few Hz, with resonant frequency of the order of 900 Hz (it 
should be observed, however, that with recent improvement in the electronics, EXPLORER 
has now a useful band of the order of a few tens of Hz). Fig. 1 shows the sensitivity curve 
of NAUTILUS. 




Figure 1: The sensitivity curve of NAUTILUS (courtesy of the NAUTILUS collaboration) 



8 



Stochastic backgrounds of GWs 




Frequency [Hz] 



Figure 2: The planned VIRGO sensitivity curve. 
4.2 Ground-based interferometers 

The first generation of large scale interferometers is presently under construction. LIGO con- 
sists of two detectors, in Hanford, Washington, and Livingston, Louisiana, with 4 km arms. 
VIRGO is under construction near Pisa, and has 3 km arms. Somewhat smaller are GEO600, 
near Hannover, with 600m arms and TAMA300 in Japan; these smaller interferometers also 
aim at developing advanced techniques needed for second-generation interferometers. Fig. ^ 
shows the planned VIRGO sensitivity curve. Interferometers are wide-band detectors, and 
will cover the region between a few Hertz up to approximately a few kHz. Comparing the 
data in fig. || with eq. (|^) we see that used as single detectors, VIRGO and LIGO could 
measure 

hln^„>10'^ . (33) 

As discussed in the previous section, much more interesting values can be obtained correlating 
two interferometers, an interferometer and a bar, or two bars. The values for one year of 
integration time are given in table 2. 



detector 1 


detector 2 




LIGO-WA 


LIGO-LA 


5 X lO^** 


VIRGO 


LIGO-LA 


4 X 10-6 


VIRGO 


LIGO-WA 


5 X 10-6 


VIRGO 


GEO 


5.6 X 10-6 


VIRGO 


TAMA 


1 X 10"^ 


VIRGO 


AURIGA 


4 X 10-^ 


VIRGO 


NAUTILUS 


7 X 10-^ 


AURIGA 


NAUTILUS 


5 X 10-'^ 



Table 2: The sensitivity for various two-detectors correlations. 



9 



4.3 Advanced interferometers 

The interferometers presently under construction are the first generation of large-scale inter- 
ferometers, and second generation interferometers, with much better sensitivities, are under 
study. In particular, the first data from the initial LIGO are expected by 2002, and the 
improvements leading from the initial LIGO to the advanced detector are expected to take 
place around 2004-2006. 

The overall improvement of LIGOII is expected to be, depending on the frequency, one 
or two orders of magnitude in hf. This is quite impressive, since two order of magnitudes in 
hf means four order of magnitudes in hlflg„{f) and therefore interesting sensitivity, h^Qg^ 
of order 10~^, even for a single detector, without correlations. 

Correlating two advanced detectors one could reach extremely interesting values: the 
sensitivity of the correlation between two advanced LIGO is estimated to be [^j 

/i^figw ^ 5 X 10"" . (34) 

4.4 The space interferometer LISA 

This project, approved by ESA and NASA, is an interferometer which will be sent into space 
around 2010. Going into space, one is not limited anymore by seismic and gravity-gradient 
noises; LISA could then explore the very low frequency domain, 10""^ Hz < / < 1 Hz. At 
the same time, there is also the possibility of a very long path length (the mirrors will be 
freely floating into the spacecrafts at distances of 5 x 10^ km from each other!), so that 
the requirements on the position measurement noise can be relaxed. The goal is to reach 
a strain sensitivity hf = 4 x 10~^^ Hz"^*^^ at / = 1 mHz. At this level, one expects first 
of all signals from galactic binary sources, extra-galactic supermassive black holes binaries 
and super-massive black hole formation. Concerning the stochastic background, going to 
low frequencies provides a terrific advantage: eq. (^ tells that 

figw(/) ~ fS,{f) . (35) 

The ability to perform an interferometric measurement is encoded into the spectral density 
of the noise, Sn{f), and with a single detector we can measure Sh{f) ~ 5'„(/). Then, if we 
are able to reach a good value of Sn at low frequency, eq. (|35| ) shows that we are able to 
measure a very small value of fig^, thanks to the factor f^. Going down by a factor 10^ 
in frequency, say from 10 Hz to 1 mHz, the factor provides an improvement by a factor 
10^^! Even used as a single detector, LISA can therefore reach very interesting sensitivities. 
Indeed, in terms of /ig^gw the sensitivity of LISA corresponds to 

hln^^if = 1 mHz) ~ 1 X 10-^2 . (36) 

We can now also understand why it is instead very difficult to build a good detector of 
stochastic GWs at high frequencies, / 3> 1 kHz. One can imagine to build an apparatus 
that measures extremely small displacements, at high frequencies, so that at some / ^ 1 
kHz, the spectral density of the noise Sn{f) is very small. However, the quantity which 
is relevant for the comparison with theoretical prediction, and on which strong theoretical 
bounds exist at all frequencies, is h^Q^^, and once one translates the sensitivity in terms of 
/igfigw, it will be very poor because the factor now works in the wrong direction. 
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5 Bounds on h^Qg^ 

In Fig. ^ we show the most relevant bounds on stochastic GW backgrounds, together with 
the experimental sensitivities discussed above. On the horizontal axis we cover a huge range 
of frequencies. The lowest value, / = 10^^*^ Hz, corresponds to a wavelength as large as the 
present Hubble radius of the universe; the highest value shown, / = 10^^ Hz, has instead 
the following meaning: if we take a graviton produced during the Planck era, with a typical 
energy of the order of the Planck or string energy scale, and we redshift it to the present 
time using the standard cosmological model, we find that today it has a frequency of order 
10^^ or 10^^ Hz. This value therefore is the maximum possible cutoff of spectra of GWs 
produced in the very early universe. The maximum cutoff for astrophysical processes is of 
course much lower, of order 10 kHz [|^ §]. So this huge frequency range encompasses all the 
GWs that can be considered. Let us now discuss in turn the various bounds. 




Figure 3: The bounds from nucleosynthesis (horizontal dashed lines, for = 4 and for 
Nu = 3.2), from COBE and from ms pulsars, together with the sensitivity of a single ground 
based detector of the VIRGO or LIGO type (blue sohd line), a bar- interferometer correlation 
(pink star), a 2-interferometer correlation (blue solid line; LIGO-LIGO, LIGO-VIRGO and 
VIRGO-GEO all give very similar sensitivities), two advanced LIGO correlation (blue dashed 
line ) and LISA (pink dashed line). 
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5.1 Nucleosynthesis bound 

The outcome of nucleosynthesis depends on a balance between the particle production rates 
and the expansion rate of the universe, measured by the Hubble parameter H. Einstein 
equation gives H"^ ~ Gjsip, where p is the total energy density, including of course Pg^. 
Nucleosynthesis successfully predicts the primordial abundances of deuterium, '^He, "^He and 
^Li assuming that the only contributions to p come from the particles of the Standard 
Model, and no GW contribution. Therefore, in order not to spoil the agreement, any further 
contribution to p at time of nucleosynthesis, including the contribution of GWs, cannot 
exceed a maximum value. The bound is usually written in terms of an effective number of 
neutrino species iVj,, and, applied to GWs, reads 

{•f=oo 

/ _ d{\og f) hln^M) < 5.6 X 10-%N, - 3) . (37) 

The limit on N^, is subject to various systematic errors, which have to do mainly with the 
issues of how much of the observed ^He abundance is of primordial origin, and of the nuclear 
processing of ^He in stars. Different bounds on N^, have therefore been proposed. The 
analysis of ref. ||^ gives Ni, < 3.2 at 95% c.l.; this is shown in fig. ^, together with a more 
conservative bound Njy < 4. Note that the nucleosynthesis bound is really a bound over the 
total energy density, i.e. over the integral of h^Qg^lf) over dlogf. However, if the integral 
cannot exceed this value, also its positive integrand cannot exceed it over a sizable range 
of frequencies. The actual bound on hlflg„{f) depends on its frequency dependence, and if 
for instance hlQg„{f) is approximatly constant between two frequencies fmm and /max, the 
bound on h^Qg^ is stronger by a factor ~ log(/max//min)- Of course, the nucleosynthesis 
bound applies only to GW produced before nucleosynthesis, and is not relevant for GW of 
astrophysical origin. 



5.2 Bounds from millisecond pulsars 

Millisecond pulsar are an extremely impressive source of high precision measurements [p!0| . 
For instance, the observations of the first msec pulsar discovered, B1937+21, after 9 yr of 
data, give a period of 1.557 806 468 819 794 5 ± 0.000 000 000 000 000 4 ms. As a consequence, 
pulsars are also a natural detector of GWs, since a GW passing between us and the pulsar 
causes a fluctuation in the time of arrival of the pulse, proportional to the GW amplitute 
hc{f). If the uncertainty in the time of arrival of the pulse is e and the total observation 
time is T, this 'detector' is sensitive to hc{f) ~ e/T, for frequencies / ~ 1/T. The highest 
sensitivities can then be reached for a continuous source, as a stochastic background, after 
one or more years of integration, and therefore for / ~ 10~^ — 10"*^ Hz. Based on the data 



from PSR B1855+09, ref. ^ gives a limit, at / = /, = 4.4 x lO^^ Hz (at 90% c.l.), 

hlng^if,) < 4.8 X 10"^ . (38) 
Since the resolution on hc{f) is proportional to 1/T and h^Qg^ ~ hi, the bound for f > f* 



IS 



hl^gM) < 4.8 X 10-9 ] , (39) 



12 



Stochastic backgrounds of GWs 



and therefore it is quite significant (better tlian the nucleosynthesis bound) even for / ~ 
10^/,= . For / < instead, the pulsar provides no limit at all. With the observation time the 
bound will improve steadily and will move toward lower and lower frequencies. It is given 
by the wedge-shaped curve in fig. ||. 



5.3 Bound from COBE 



Another important constraint comes from the COBE measurement of the fluctuation of the 
temperature of the cosmic microwave background radiation (CMBR). The basic idea is that 
a strong background of GWs at very long wavelengths produces a stochastic redshift on 
the frequencies of the photons of the 2.7K radiation, and therefore a fluctuation in their 
temperature (Sachs- Wolfe effect). The analysis of refs. [jl2|, |^, where the effect of multipoles 
with 2 < / < 30 is included, gives a bound 



hln^if) < 7 X 10- 




(3 X 10'^^ Hz < / < 10^^^ Hz). 



(40) 



This bound is stronger at the upper edge of its range of validity, / ~ SOHq ~ 10 Hz, 
where it gives 

hln^if) < 7 X 10-1^ , (/ ~ 10-'' Hz) . (41) 



6 Theoretical predictions 

Many cosmological production mechanisms have been proposed in recent years. Four exam- 
ples are shown in fig. ^ 

Inflation. The amplification of vacuum fluctuations at the transition between an in- 
flationary phase and the radiation dominated era produces GWs shown as the curve (a) in 
fig. ^ and is one of the oldest |TB] and most studied examples. The condition that the COBE 
bound is not exceeded puts a limit on the value of h^Qg^ at all frequencies, which is below 
the experimental sensitivities, even for LISA. So, while it is one of the best studied examples, 
it appears that it is not very promising from the point of view of detection. 

String cosmology. In a cosmological model which follows from the low energy action 



of string theory [|1J, |15] the amplification of vacuum fluctuations can give a much stronger 
signal ||16| , p!7| . The model has two free parameters that reflect our ignorance of the large 
curvature phase. The curve (b) in fig. ^ shows the very interesting signal that could be 
obtained for some choice of these parameters, while for other choices the value of h^Qg^ at 
VIRGO/LIGO or at LISA frequencies becomes unobservably small. 



ries 



Cosmic strings. These are topological defects that can exist in grand unified theo- 
TB|, and vibrating, they produce a large amount of GWs, shown in curve (c) of fig. H. 



Cosmic strings are characterized by a mass per unit length /i, and the most stringent bound 
on GAT/i comes from msec pulsars, and it is of order 10~^. 



Phase transitions. Another possible source of GWs is given by phase transitions in 
the early universe. In particular, a phase transition at the electroweak scale would give a 
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signal just in the LISA frequency window, while the QCD phase transition is expected to 
give a signal peaked around / = 4 x 10~^ Hz. However, the signal is sizable only if the phase 
transition is first order and, unfortunately, in the Standard Model with the existing bounds 
on the Higgs mass, there is not even a phase transition but rather a smooth crossover, so 
that basically no GW is produced. However, in supersymmetric extensions of the Standard 
Model, the transition can be first order, and a stronger signal could be obtained. Depending 
on the strength of the transition, one could even get a signal such as curve {d) of fig. |. 




-18.0 -14.0 -10.0 -6.0 -2.0 2.0 6.0 10.0 
Log[f(Hz)] 

Figure 4: The backgrounds predicted, with optimistic choice of parameters, by (a) inflation, 
(b) string cosmology, (c) cosmic strings, (d) a first order phase transition at the electroweak 
scale, together with the bounds and sensitivities of fig. |. 



Finally, there are very interesting astrophysical backgrounds, coming from a large number 
of unresolved sources. These are displayed in fig. |[ For a discussion, see |]19[ and the contri- 
bution of Raffaella Schneider to these proceedings. Another important issue, especially for 
LISA, is also how to discriminate cosmological from astrophysical backgrounds, see eg. . 



The conclusion that emerges looking at fig. |^ is that in the next few years, with the first 
generation of ground based interferometers, we will have the possibility to explore five new 
order of magnitude in energy densities, probing the content in GWs of the universe down 
to /ip^gw ~ 10~^. At this level, the nucleosynthesis bound suggest that the possibility of 
detection are quite marginal. It should not be forgotten, however, that nucleosynthesis is a 
(beautiful) theory, with a lot of theoretical input from nuclear reaction in stars, etc., and 
its prediction is by no means a substitute for a measurement of GWs. With the second 
generation of ground based interferometers and with LISA, we will then penetrate quite 
deeply into a region which experimentally is totally unexplored, and where a number of 
explicit examples (although subject to large theoretical uncertainties) suggest that a positive 
result can be found. 
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Stochastic backgrounds of GWs 




10-6 lQ-5 lQ-4 lQ-3 0.01 0.1 1 10 102 IQQ IQ4 

Frequency [Hz] 



Figure 5: Stochastic backgrounds of astrophysical origin. From left to right, black hole- 
black hole binaries, white dwarf-white dwarf binaries, neutron star-neeutron star binaries 
(the curve extending up to 10^), a pregalactic Population III stars, r- modes of neutron stars 
and supernova collapse to black holes (figure kindly provided by Raffaella Schneider). 
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